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INTRODUCTION 


The  elastic-plastic  stress-strain  relations  for  a  rod  have  been  derived 
by  T.  Wright  (ref  1),  and  the  differential  equations  of  the  rod  itself  have 
been  cast  into  variational  form.  The  variational  principle  using  bilinear 
functional  and  adjoint  formulation  has  served  as  a  basis  to  determine 
numerical  solutions  by  the  finite  element  method.  This  principle  has  now  been 
extended  to  coupling  systems  such  as  in  the  impact  dynamics.  The  present 
hyperbolic  type  partial  differential  equation  has  two  dependent  and  two 
independent  variables,  with  coupling  in  the  spatial  domain.  This  new 
formulation  is  also  ready  to  be  used  for  the  coupled  impact  problem  which  is 
given  in  the  Appendix  of  this  report. 

THE  VARIATIONAL  PRINCIPLE 

A  dynamical  system  can  be  modeled  by  the  matrix  vector  partial 
differential  equation. 

L(0  y(0  -  -Q(  O  (1) 

with  appropriate  boundary  and  initial  conditions.  In  the  above  equation,  L  is 
a  matrix  linear  operator  in  both  spatial  and  temporal  domain,  y  is  a  vector 
dependent  variable,  Q  is  a  vector  forcing  function,  and  £  represents  all 
independent  variables,  both  spatial  and  temporal. 

The  inner  product  <  >  of  an  adjoint  forcing  function  Q  and  the  solution 
(y(O)  of  Eq.  (1)  can  be  used  for  the  purpose  of  estimation.  This  inner 
product  is  <Q,y>. 

^■Wright,  Thomas,  W.,  "Nonlinear  Waves  in  Rods,"  Technical  Report 
ARBRL-TR-02324 ,  U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen,  MD  21005, 
May  1981. 
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An  accurate  estimation  can  be  made  by  constructing  a  variational 
principle  (refs  2-5).  By  using  the  adjoint  variable  y  as  a  Lagrange 
multiplier  for  Eq.  (1)  and  adding  to  <Q,y>,  we  have 

Ji[y>y]  =  <Q,y>  +  <y,(Q+Ly)>  =  <Q,y>  +  <y,Q>  +  <y,Ly>  (2) 

To  keep  the  system  symmetrical,  let  us  define  the  adjoint  system  as 

L(Qy(0=-Q(Q  (3) 

By  using  the  original  variable  y  as  a  Lagrange  multiplier  for  Eq.  (3)  and 
adding  to  <Q,y>,  we  have 

j2[y,y]  =  <Q»y>  +  <y»(Q+Ly)>  =  <Q»y>  +  <y>Q>  +  <y»Ly>  (4) 

The  relationship  of  the  adjoint  system  to  the  original  system  is 

A  - 

*  D  =  <y,Ly>  -  <y,Ly>  =  De  (5) 

where  D  is  the  bilinear  concomitant  (ref  2).  Combining  Eqs.  (2),  (4),  and  (5) 
one  obtains 

JX  =  J2  +  De  (6a) 

In  order  to  keep  the  functional  symmetrical,  we  have 

A  1 

Jl  =  -  [Ji  +  ^2  +  Del  (6b) 


^Stacey,  Weston,  M.,  Jr.,  Variational  Methods  in  Nuclear  Reactor  Physics, 
Academic  Press,  1974. 

^Shen,  C.  N.  and  Wu,  Julian  J.,  MA  New  Variational  Method  for  Initial  Value 
Problems,  Using  Piecewise  Hermite  Polynomial  Spline  Functions,”  ARO  Report 
81-3,  Proceedings  of  the  1981  Army  Numerical  Analysis  and  Computers 
Conference,  1981. 

^Shen,  C.  N. ,  "Variational  Principle  for  Gun  Dynamics  With  Adjoint  Variable 
Formulation,"  Proceedings  of  the  Third  US  Army  Symposium  on  Gun  Dynamics, 
Volume  II,  May  1982,  p.  IV-108. 

^Shen,  C.  N. ,  "On  the  Extremum  of  Bilinear  Functional  for  Hyperbolic  Type 
Partial  Differential  Equations,"  ARDC  Technical  Report  ARLCB-TR-84008 ,  Benet 
Weapons  Laboratory,  Watervliet,  NY,  April  1984. 
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which  is  of  the  form 


Similarly , 


Ji  =  <Q,y>  +  <y,Q>  +  -  <y,Ly>  +  -  <y.Ly>  +  " 

1  -  1  —  De 

j2  =  <Q,y>  +  <y,Q>  +  -  <y.Ly>  +  ~  <y»Ly>  — ~2 


(6c) 


(6d) 


INTEGRAL  OF  BILINEAR  EXPRESSION 

The  integral  of  a  bilinear  expression  for  a  two-dimensional  problem 
having  a  system  of  second  order  partial  derivatives  in  time  and  in  space  can 
be  written  as 

xb  tb 

I=f  J  fl[y(x,t),y(x,t)]dtdx  (7) 

xo  to 

where  J2[y,y]  is  a  given  bilinear  expression  in  the  form 

fi[y,y]  =  ytTgyt  "  yxlAyx  -  yTpy  +  yxTry  +  yTyyx  (8) 

The  subscripts  t  and  x  indicate  the  partial  derivatives  for  the  functions  y 
and  y.  The  matrices  A,  B,  and  P  are  diagonal  and  T  and  N  are  off-diagonal. 
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Equation  (8)  can  be  integrated  by  parts.  Two  different  forms  of 
integration  and  end  conditions  are  given.  The  first  form  of  the  integral  is 
obtained  by  integrating  by  parts  on  the  adjoint  variable. 


3 


(11) 


which  gives 


la  =  -  //{yTBytt  -  yTAyxx  +  yTpy  +  yTpyx  -  yxNyxMtdx 

+  /  yTBytdx  -  /  yXAyxxdt  +  /  yTPydt 
Ia  =  If  -  w(b2iwtt-ai  ^xx+Pl  2w+Ylux-niux)dtdx 

+  //  -  u(b22utt_a22uxx+P22u+^2w2-n2wx)dxdt 

xb  -  ,  ~ 

+  J  [wb^^Wj-  +  ubnij;]  dx  + 
xo  to 

xb  -  xb  tb  -  xb 

/  [wa^x  +  032^x1  dt  +  J  [WY]U  +  uY2w]  dt 


(12) 


xb  tb 

I  =  -<y,Ly>  +  J  [bi^rot  +  b2^iut]  dx 


tb  -  -  xb  tb  -  -  xb 

-  /  [ai  ^wwx  +  a2  niiix]  dt  +  J  [Yiwu  +  Y2uw]  dt 


(13) 


On  the  other  hand,  we  can  perform  integration  on  the  original  variable  to  give 

ib  =  -J/{yXBytt  “  yTAyxx  +  yXpy  -  yXiqTx  +  yTNTyx}dtdx 


/yxBytdx  -  /ypAyxdt  +  / y^N^yd  t 


(14) 


which  gives 


Ib  =  //(■*w)(bi2wtt-  al  ^xx+Pl  2W“  Y2UX+  h2ux)  d  t  dx 

+  11  (-u)(b12utt-a22uxx+P22u-Ylwx+niwx)dtdx 

xb  -  -  fcb  tb  -  -  xb  tb  _  -  xb 

+  /  [wbi2wt+ub22ut]  dx  -  f  [wa^  ^fuamj  dt  +  J  [wn2u  +  uniw]  dt  (15) 


Ib  =  -<y ,Ly>  +  /  [bi2wwl-  +  b22uut;]  dx  -  f  [ al +  a22uuxl  dt  + 
xo  tQ  tQ  x0 

tb  -  -  xb 

+  |  [  r^wu  +  ri]_uw]  dt 


-  tb 


tb 


-  xb 


(16) 
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To  keep  the  form  symmetrical,  we  take  the  average  of  the  previous  two 


expressions 


1  1  xb  tb  1  ~ 

I  =  -  Ia  +  -  Ib  =  -J  J  -  (yLy+yLy)dtdx  + 

2  a  2  b  xc  t0  2 


+  J(yTByt+yTByt) I  bdx  -  -  J(yTAyx+yTAyx) I  bdt  +  -  J(yTPy+yTNTy)dt  (17) 


to  2 


x0  2 


which  gives 


1  -  1  --  1  xb  ,  -  -  9 

I  - - <y  ,Ly> - <y,Ly>  +  -  J  [bi  z(wwt+wwt)  +  b2  nuut+uut)  ]  dx 

2  2  2  x0  tQ 

1  tb  -  -  -  -  xb 

-  -  J  [a^ z(wwx+wwx)  +  b2  ^(uujj+uujj)  ]  dt 

2  t0  x0 

1  tb  -  -  xb 

+  -  J  [(Yl+ni)wu  +  (Y2+n2)uwl  dt 

2  tn  x_ 


(18) 


where 


_ 

!  LU 

1*12  I 

-  j  mi 

1*12  1 

L  =  | 

1 

,  L  =  1  _ 

—  ! 

1  1*1 

l22  j 

1  1*21 

1 

l22  1 

(19) 


L11  =  L11 


32  32 

H2  +  pi2  -  ai2  ”2 


32 

,  L22  =  L22  =  b22  ---  +  P22 

dt * 


l12  =  l12 


3 

3x 


l21  =  l21  “  ” 
3x 


32 

a22  — 2  (20) 


(21) 


It  is  noted  that  we  have  used  the  following  relationship  in  obtaining  Eq. 

(21) 

hi  -  Yl  =  1  and  h2  ~  Y2  =  “1  (22) 

by  comparing  Eqs.  (10),  (21),  and  the  last  equation  in  the  Appendix. 

For  a  spatial  and  temporal  partial  system,  Eq.  (5)  becomes 

xb  tb  —  xb  tb  — 

D  =  /  /  yLydtdx  -  /  /  yLydtdx  (23) 

X°  tQ  xQ  t0 
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By  equating  Eqs.  (11)  and  (14)  and  solving  for  D  in  Eq.  (23),  we  are 
converting  the  double  integral  into  single  integrals  in  terms  of  the  boundary 


conditions. 

We  can  express  the  quantity  De  as  the  sum  of  three  parts  for  end 
conditions  D2,  and  D3  as 

De  =  +  D2  +  D3 

The  terms  in  involve  the  initial  conditions  of  y  and  y  as 

xb  “  tb 

Dl  =  /  (yTByt-yTByt) I  dx 


(24) 


xb  „  -  -  2  -  tb 

=  J  [bi  (wwt-wwt)  +  b2  (uut-uut)]  dx 


(25) 


Ao  o 

The  terms  in  D2  involve  the  boundary  conditions  from  the  second  partials 
of  y  and  y  as 

tb  -  xb 

D2  =  -/  (yTA.yx-yTAyx)  I  dt 


ftb  2  -  -  -  xb 

-  J  [ai  (wwx-wwx)  +  a2z(uux-uux)]  dt 

Xr\ 


(26) 


The  terms  in  D3  involve  the  boundary  conditions  from  the  coupling  terms. 

/tb  —  -  .cb  -  -  xb 

D3  =  J  (yTTy-yTNTy)dt  =  J  [(Yl~ni)wu  +  (Y2-h2)uw]  dt 


(27) 


THE  SYMMETRICAL  ADJOINT  SYSTEM 


The  adjoint  independent  variable  x  in  Figure  1  can  be  expressed  as 

Tb  -  T  t  -  t0 
Tb  “  To  bb  —  tQ 


(28) 


6 


**“>  c?> 


X 


-h)  ~  Lc 


utx,U£b)= 

(^t-%)-  ~^ttx/  i  =&>) 


Figure  1.  Image  Reflection  of  the  Adjoint  System. 


which  gives 


and 


T  =  Tb  for  t  =  tQ 
T  =  T0  for  t  =  tb 
It  is  noted  from  Eq.  (28)  that 

Tb  —  To  -  Tb  —  *-o 

T  "  Tb  +  co  "  c 
d  T  =  -dt 

d  _  d 
dt  dt 

and 


(29a) 

(29b) 

(30a) 

(30b) 

(30c) 

(30d) 

(30e) 


y(x,t)  =  y(x,x  =  xb+t0-t) 

Let  us  assume  that  the  adjoint  system  shown  in  Figure  1  gives 

A  A 

y(x,t=t)  =  y(x,t=tb+tQ-t)  (31a) 

A  A 

yt(x,t=t)  =  -yt(x,t=tb+t0-t)  (31b) 

A 

yx(x,t=t)  =  yx(x,t=tb+tQ-t)  (31c) 

where  t  is  a  dummy  variable  for  t. 

We  may  define  the  adjoint  system  as  the  image  reflection  in  the  time 
domain  of  the  original  system.  Equation  (31)  yields  the  following  known 
initial  conditions: 

A  A 

y(x,t=tb)  =  y(x,t=tQ)  (known)  (32a) 

A  A 

yt(x,t=tb)  =  -yt(x,t=t0)  (known)  (32b) 

The  interpretation  of  the  above  equations  gives  the  initial  conditions  of  the 
original  system  as  the  far  end  conditions  for  the  adjoint  system,  since  the 
adjoint  system  is  a  reflected  mirror  of  the  original  system  in  time. 
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INITIAL  CONDITIONS  FOR  THE  ADJOINT  SYSTEM 

We  take  a  symmetrical  approach  for  the  initial  conditions  of  the  adjoint 

system  as 

y(x,t=tb)  =  y(x,t-t0)  ,  yt(x,t-tb)  =  -yt(x>t=to)  (33) 

y(x,t=t0)  =  y(x,t-tb)  ,  yt(x,t-t0)  =  -yt(x,t-tb)  (34) 

where  y  and  y  are  given  in  Eq.  (10).  Thus  Eq.  (25)  becomes 
xb 

D1  =  J  b12dx{[w(x,t=t0)wt(x,t=tb)  +  w(x,t=tb)wt(x,t=t0) ] 
xo 

-  [w(x,t=tb)wt(x,t=t0)  +  w(x,t=t0)wt(x,t=tb)] } 
xb 

+  /  b22dx([u(x,t=t0)ut(x,t=tb)  +  u(x,t=tb)ut(x,t=t0)] 
xo 

-  [u(x,t=tb)ut(x,t=t0)  +  u(x,t=t0)ut(x,t=tb)] }  =  0  (35) 

Since  the  integrand  of  Eq.  (35)  is  zero,  the  above  satisfies  Eq.  (25). 

THE  GENERALIZED  BOUNDARY  CONDITIONS 

Let  us  consider  the  operator  L  in  Eqs.  (19)  through  (21) •  It  is  assumed 
that  elastic  springs  are  installed  at  the  ends  such  that 

yx(xb»t)  =  Kby(xb,t)  ,  yx(xb,t)  =  Kby(xb,t)  (36a) 

yx(x0,t)  =  -K0y(x0,t)  ,  yx(xQ,t)  =  -K0y(x0,t)  (36b) 

where  Kb,  Kq  are  diagonal  matrices.  If  Eq.  (36)  is  substituted  into  Eq.  (26), 
we  have 

D2  =  0  (37) 

Since  D^  =  D2  =  0,  Eq.  (24)  becomes  D3  as  given  in  Eq.  (27) 

De  =  D3  (38) 
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CONDITIONS  FOR  THE  COUPLING  TERMS 

The  sum  of  the  functionals  Jb  +  I  is  obtained  by  adding  Eqs.  (6c)  and 
(18)  as 

Xb  tfc  _  1 

Jl  +  l  +  f  j  (Qy+yQ)dxdt  +T+B+V+-D3 


where 


xo  to 


and 


From  the  last  two  equations,  one  obtains 


1  tb  xb 

V  H —  D3  =  J  [  (  Yiuw+Y2wu)  ]  dt 


(39) 


1  xb 

[t>i2 

-  -  tb 

(wtw+wtw)  +  b2  (utu+utu)]  dx 

(40) 

2  xo 

-  ;  Jtb 

tai 

to 

-  -  xb 

2(wxw+wxw)  +  a 2  z(uxu+uxu) ]  dt 

(41) 

2  tQ 

1 

tb 

xo 

-  xb 

(42) 

V  -  - 

J 

[  (  Yl+hi)uw  +  (  Y2Tl2)wul  dt 

2 

xQ 

1 

tb 

-  xb 

(43) 

d3  -  r 

I 

[(Yi-hi)uw  +  (Y2-D2)vm]  dt 

2 

xo 

(44) 


We  can  let  Eq.  (44)  vanish  by  choosing 

Yi  =  Y2  =  0 

which  gives 


From  Eq.  (22)  one  obtains 


V  +  -  Do  =  0 
2  J 


hi  =  1  +  Yl  =  1 
T12  =  1  +  Y2  *  “1 


(45) 

(46) 

(46a) 

(46b) 
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Thus,  the  functional  for  the  original  variables  and  adjoint  variations  becomes 

Jl  =  -I  +  <Q,y>  +  <Q,y>  +  T  +  B  (47) 


The  sum  of  the  two  functionals  J2  +  I  is  obtained  by  adding  Eqs.  (6d)  and  (18) 
as 


.xb  tb  -  -  1 

J 2+I=f  f  (Qy+yQ)dxdt  +T+B+V--D3 

xo  to  2 


where  T,  B,  V,  and  (1/2)D3  are  given  in  Eqs.  (40)  through  (43).  By 

subtracting  (1/2)D3  from  V  we  have 

1  tb  -  -  xb 

V  -  -  D3  =  /  [ ( niuw+n2wu) ]  dt 

2  tQ  xQ 


(48) 


(49) 


In  this  case  we  let 

ni  =  q2  =  0  (50) 

Then  from  Eq.  (22)  one  obtains 

yi  -  -i  +  m  =  -i  (5i) 

Y2  =  1  +  n2  =  1  (52) 

Thus,  the  functional  for  the  adjoint  variables  and  original  variations  becomes 

J2  -  -I  +  <Q,y>  +  <Q,y>  +  T  +  W  (53) 

which  gives  the  same  form  as  shown  in  Eq.  (47). 


THE  FIRST  VARIATION 


By  taking  the  variations  <5y  and  6y  separately,  we  let 

SJ  =  6J (  6y)  +  5J2(  6y)  =  0  +  0  (54) 

Then  one  obtains  from  Eqs.  (40),  (41),  and  (47)  that 

SJi(Sy)  =  -6l(6y)  +  //  Q6y  dxdt  +  5T(  6y)  +  <$B(  6y)  =  0  (55) 

where 

-  1  xb  -  t^ 

ST(Sy)  -  -  /  [bi  2(wt  6w+w6wt)  +  b2  2(ut  6u+uSut)  ]  dx  (56) 

2  xo  tQ 
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(57) 


1  tb  -  -  -  -  xb 

6B(6y)  =  -  -  /  [aq2(wx6w+w6wx)  +  a2  2(ux<Su+u6ux)  ]  dt 

2  t0  xo 

and  -6l(6y)  can  be  derived  from  Eq.  (18)  with  Yq  =  T2  =  0  and  nl  = 

-  Xb  tb  -  -  2 

-6l(6y)  =  -/  J  {[bi^tfiw,.  -  pi^dw  -  aq  wx<5wx]  + 
xo  to 

[b22ut6ut  -  P22u6u  -  a22ux6ux]  +  [ux<Sw-wx6u]  }dtdx 
The  second  term  on  the  right  side  of  Eq.  (54)  is 

6J2( 6y)  =  -6I(6y)  +  //  Q6y  dxdt  +  <ST(6y)  +  <$B(6y)  =  0 


h2 


1 


(58) 


(59) 


where 


1  xb  -  -  -  tb 

6l(  6y)  =  -  /  [bq2(w6wt+wt  <Sw)  +  b22(u6ut:+ut:6u)  ]  dx 

2  x„ 


(60) 


—  —  —  —  x  b 

6B (6y)  =  -  -  J  [aq2(w6wx+wx6w)  +  a2  2(u6ux+ux6u)  ]  dt 
2  tn  xo 


(61) 


and  6I(6y)  can  be  derived  from  Eq.  (18)  with  hq  -  h2  -  0  and  — Yl  -  Y2 


-6l(6y)  =  -  J  J  {[bq^j-dwj-  -  pq^dw  -  aq^^w^  + 
xo  co 

[b22U(-6uj.  -  p22u6u  -  a22ux6ux]  +  [— wx6u  +  ux6w]  }dtdx  (62) 

The  adjoint  equation  has  the  same  form  of  the  original  equation  by  dropping 
and  adding  the  bars  simultaneously  on  every  variable. 

Equations  (55)  through  (58)  are  the  key  equations  to  be  used  for  the 
finite  element  method.  It  is  noted  that  the  first  variation  &Jq( 6y)  is  the 
'same  as  the  first  variation  6J2( 6y)  by  adding  or  dropping  the  bar  on  top  of 
the  variables  and  their  variations.  We  do  not  need  to  solve  for  the  adjoint 
system  in  Eqs.  (39)  through  (41)  since  these  give  exactly  the  same  solutions 
as  the  ones  of  the  original  system. 
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CONCLUSIONS 


The  functional  in  bilinear  matrix  vector  form  is  symmetrical  about  the 
original  variables  and  the  adjoint  variables.  The  Euler- Lagrange  equations 
for  the  coupling  systems  are  derived  using  the  fundamental  lemma  of  the 
calculus  of  variations.  By  integrating  the  bilinear  matrix  vector  expression 
by  parts,  one  can  obtain  the  bilinear  concomitant  in  terms  of  initial  and 
boundary  terms.  The  adjoint  system  can  be  arranged  in  a  manner  that  it  is  a 
reflected  mirror  of  the  original  system  in  time.  Thus  the  initial  conditions 
for  the  bilinear  concomitant  become  zero.  Generalized  boundary  conditions 
using  many  types  of  "springs"  relating  the  various  spatial  partial  derivatives 
can  be  defined  to  satisfy  the  boundaries  of  the  concomitant.  Algorithms  are 
developed  for  use  in  the  finite  element  method  by  taking  the  first  variations 
of  the  functional.  These  algorithms  are  simplified  because  the  adjoint  system 
gives  exactly  the  same  solutions  as  those  of  the  original  system. 
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APPENDIX 


The  wave  equation  in  rods  derived  by  T.  W.  Wright  (ref  A-l)  is  given  as 
the  following  system: 

2  X  c2 

(A-la) 


2  X 

+  77;:  u5 


X+2  y 


c 

Cl 


-  2  WTT 


X+y  X  c2 

UPP  -  [5 - u  +  4  -  w^]  =  ---  uxx 

y  y  c2z 


which  can  be  transformed  to 


X+2  y  X+2  y 

wtt  -  —  wxx)  -  UX  =  0 


(A-lb) 


(A-2a) 


y  c2  y  2(£+y) 

wx  +  (77  2  utt  ”  77  uxx  +  i  u)  -  0 

4X  c2z 


4  X 


(A-2b) 


With  appropriate  group  of  parameters,  we  have  the  following  form: 

^l^tt  "  al^wxx  +  PI  ^  "  ux  =  0  (A-3a) 

wx  +  b22utt  -  a22uxx  +  p22u  =  0  ( A-3b) 

The  above  system  of  equations  can  be  expressed  by  a  matrix  vector  form  of 
equations  as 


_i  r  -j 

L11  l12  I  j  w  j 

l21  l22  I  I  u  I 

I  I  I 


'“I 
0  I 
I 

0  I 


(A-4) 


where 


a2  a2  a 

Ml  -  M  --2  -  n  -^2  +  Pi  -  M2  -  -  - 


(A-5) 


^"^-Wright,  Thomas  W.,  "Nonlinear  Waves  in  Rods,”  Technical  Report 

ARBRL-TR-02324,  U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen,  HD 
21005,  May  1981. 
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(A-6) 


L21  =  3 x  and  L22 


b22  ~2  -  a2 


2  --2  + 

3x2 


P2‘ 


which  can  be  written  as  Eq.  (1)  in  the  text. 

The  notations  for  the  wave  equation  in  rods  are: 
w  =  axial  displacement 
u  =  radial  strain 

5  =  z/a  =  nondimens ional  axial  coordinates 
t  =  ct/a  =  nondimensional  time 
a  =  elastic  stored  energy  per  unit  length 
c^  =  /(X+2p)/p  =  longitudinal  wave  speed 
C2  =  / h/P  =  shear  wave  speed 
X  and  p  are  Lame'  constants 

The  above  system  of  equations  was  first  developed  in  different  form  by 
Mindlin  and  Herrmann  (ref  A-2)  and  can  be  grouped  into  a  single  equation  as: 


32  c2  32  32  c2  32 

^3£2  c^2  3t2  3£2  C2^  3t2 


Xf  p  cb  3  2 

6 - -  (— 

p  ci2  3C2 


— 3)](w  or  u)  =  0 
cb  cb 

(A-7) 


A-2Mindlin,  R.  D.  and  Herrmann,  G. ,  "A  One-Dimensional  Theory  of  Compressional 
Waves  in  an  Elastic  Rod,”  Proceedings  of  the  First  U.S.  National  Congress 
of  Applied  Mechanics,  1950,  pp.  187-191. 
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